Maximal height statistics for 1//" signals 
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Numerical and analytical results are presented for the maximal relative height distribution of 
stationary periodic Gaussian signals (one dimensional interfaces) displaying a 1//" power spectrum. 
For < a < 1 (regime of decaying correlations), we observe that the mathematically established 
limiting distribution (Fisher-Tippett-Gumbel distribution) is approached extremely slowly as the 
sample size increases. The convergence is rapid for q > 1 (regime of strong correlations) and a highly 
accurate picture gallery of distribution functions can be constructed numerically. Analytical results 
can be obtained in the limit a oo and, for large a, by perturbation expansion. Furthermore, using 
path integral techniques we derive a trace formula for the distribution function, valid for a — 2n even 
integer. From the latter we extract the small argument asymptote of the distribution function whose 
analytic continuation to arbitrary a > 1 is found to be in agreement with simulations. Comparison 
of the extreme and roughness statistics of the interfaces reveals similarities in both the small and 
large argument asymptotes of the distribution functions. 
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I. INTRODUCTION 

Whereas the extreme value statistics (EVS) of indepen- 
dent and identically distributed (i.i.d.) random variables 
has been thoroughly understood for a long time 1, 2, _J|, 
our knowledge about the EVS of correlated variables is 
less general. Many natural processes, like flood-water 
levels, meteorological parameters, and earthquake mag- 
nitudes IE 1^1 are, however, characterized by large 
variations, a phenomenon connected to long term cor- 
relations. Since extremal occurrences in physical quanti- 
ties may be of great significance, it is essential to develop 
an understanding of EVS in the presence of correlations. 
The last few years have seen increased activity in this 
direction, with several particular cases worked out in de- 
tail. For example, extremal height fluctuations in 1-1-1 
dimensional Edwards- Wilkinson surfaces have been in- 
vestigated recently 0, 0, and a nontrivial distribution 
function, the Airy distribution, was found analytically 
for the stationary surface. Equivalently, considering the 
latter as a time signal, this result relates to maximal dis- 
placements in Brownian random walks. Other studies of 
surface fluctuations also demonstrate the effect of correla- 
tions on EVS, and several examples show that nontrivial 
EVS may emerge even in the simplest surface evolution 
models [3, ^3 lilt IT^ IT^ . Remarkable connections 
have also been found between EVS and propagating front 
solutions, exploited in such problems as random fragmen- 
tation ^3 , or random binary-tree searches . Correla- 
tions have also been shown to play an important role in 
effecting extreme events in weather records 0| . To 
summarize, problems related to extremes regularly arise. 



and it is a fundamental question whether they obey a 
limit distribution characterizing i.i.d. variables, or some 
special, nontrivial, statistics emerges. 

In order to develop an intuition about the effect of 
correlations, we shall consider here the EVS of periodic 
signals displaying Gaussian fluctuations with 1//" power 
spectra. While we shall use the terminology of time sig- 
nals, one-dimensional stationary interfaces may equally 
be imagined, with the same spatial spectrum, and peri- 
odic boundary conditions. Systems with 1//" type fluc- 
tuations are abundant in nature, with examples ranging 
from voltage fluctuations in resisto rs Jl9l| . through tem- 
perature fluctuations in the oceans |2(]| | , to climatological 
temperature records , to the number of stocks traded 
daily [22|. In addition, most of these fluctuations ap- 
pear to be Gaussian, thus our results may have relevance 
in answering questions about the probability of extreme 
events therein. 

The 1//" signals we consider are rather simple in the 
sense that they decompose into independent modes in 
Fourier space. The modes are not identically distributed, 
however, giving rise to temporal correlations, which are 
by now well understood (see Sec. |nj. Gorrelations are 
tuned by a, yielding signals with no correlations (a = 
0), decaying (0 < a < 1), and diverging correlations 
(1 < a < 00). Thus 1//" processes are also well suited 
for studying the effect of a wide range of correlations on 
extreme events in signals. 

The central quantity we investigate is the maximum 
relative height (MRH), first studied in |2l|. This is the 
highest peak of a signal over a given time interval T, 
measured from the average level. Specifically, for each 
realization of the signal, h{t), the MRH is 



hm '■— max/i(t) — h{t) . 



(1) 
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where maxt h{t) is the peak of the signal and h{t) is its 
time average. The MRH, hm, varies from realization 
to realization, and is therefore a random variable whose 
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probability density function (PDF), denoted by P{hm), 
we would like to determine. The physical significance of 
hm is obvious. For instance, in a corroding surface it 
gives the maximal depth of damage or, in general, it is 
the maximal peak of a surface. To name another exam- 
ple, when natural water level fluctuations are considered, 
it is related to the necessary dam height. 

Since the Fourier components of the signal are indepen- 
dent variables, it is relatively easy to generate hm-s nu- 
merically and thereby obtain sufficient statistics for sam- 
pling P{hra) (see Sec lIIIjl . Scanning through < a < oo 
reveals that ac — 1 separates two regions with distinct 
behaviors in both the limiting functions and the conver- 
gence to them as the signal length (0 < i < T) tends to 
infinity. At a = 0, the signal is made up of i.i.d. vari- 
ables and the EVS is governed by the FTG distribution, 
which is one of the three possible limit distributions for 
i.i.d. variables in the traditional categorization pS. l25l |. 
In fact, this property extends to the whole < a < 1 
interval |2^, where the correlations decay in a power- 
law fashion. Our results indicate that, at least in the 
< a < 0.5 region, not only the limit distribution but 
the convergence to it follows closely the logarithmically 
slow convergence which characterizes a = (Sec. IIV|I . 
We find that the convergence further slows down in the 
< a < 0.5 region and it remains an open question 
whether it is slower or not than logarithmic. 

For a > 1, the signal becomes rough, that is, the cor- 
relations diverge with signal length, and we find that the 
qualitative features of the EVS in this range are the same 
as in the a = 2 case (Sec.0), exactly solved by Comtet 
and Majumdar Namely, the divergent scale of the 

extreme values {hm) , where /3 = (a — l)/2, is pro- 
portional to the scale of the fluctuations in the signal 
(square root of the roughness in interface language) and, 
furthermore, the large- and small-argument asymptotes 
of the limiting distribution functions are of similar type. 
In order to demonstrate these similarities, we study the 
generalized, higher order, random acceleration problem 
[a = even integer) in Sec. I VII and calculate the propaga- 
tor of this process. Using this result, we develop a gener- 
alization of the trace formula fSec lVIlT) which was instru- 
mental in solving the a — 2 problem. It turns out that 
the trace formula can be written in a scaling form, which 
yields the scale of the MRH values fSec lVin)) as well as, 
under a rather mild and natural assumption, the small- 
argument asymptote of the MRH distribution fSec. IIX|I . 
Our numerical evaluations of the distributions are all in 
excellent agreement with the analytical results. 

Analytical results can also be obtained in the a — > oo 
limit (Sec. 0), where the lowest frequency mode deter- 
mines the shape of the signal. We flnd that the MRH 
distribution has the functional form ~ xexp {—x^). Cor- 
rections to the a — > oo limit may be obtained by keeping 
the lowest frequency modes. With only three modes, a 
satisfactory description of the whole a > 6 region can be 
obtained (Sec. O. Since both the a = 2 and a = oo 
results suggest that the large-argument tail of the distri- 



bution takes the form ~ x"*" exp (— a:^), we checked this 
property for other a-s as well, and found it to be an ex- 
cellent description for all a > 1. 

The common scaling properties of the maximal height 
and the root mean square height for a > \ lead us to 
compare the MRH distributions to the roughness distri- 
butions of 1//" interfaces We find in Sec. Ixj that, 
in addition to the general shape of the PDF-s, both the 
small and large argument asymptotes of these functions 
have analogous functional forms provided the replace- 
ment hm ^ {roughnessY^^ is made. Similar conclusions 
can also be reached when P{hm) is compared with the 
distribution of maximal intensities . 

Concluding remarks are collected in Sec. IXIII while de- 
tails of the calculations of the generalized random accel- 
eration process and of the large a expansion are given in 
Appendix ^ and respectively. 

II. GAUSSIAN PERIODIC 1//° SIGNALS 

We consider Gaussian periodic signals h{t) = h{t + T) 
of length T. The probability density functional of h{t) is 
given by 

7'[/i(i)] -e-^l''(*)l, (2) 

where the effective action S can be formally defined in 
real space but, in practice, is defined through its Fourier 
representation 

N/2 

^[c„;a] = 2Ari-"^n"|c„p. (3) 

n=l 

Here A is a stiffness parameter which is set to (27r)"/2 
hereafter (for the details and notation we follow |23|), 
and the c„-s are the Fourier coefficients of h{t) 

N/2 

h{t)^ ^ c„e^— (4) 

n=-N/2+l 

where c* = c_„ and their phases (for n ^ -^72) are 
independent random variables uniformly distributed in 
the interval [0,27r], while C7V/2 is real. Since cq does not 
appear in the action (jS)) we can set the average of the 
signal to zero, i.e. cq — 0. Note that the cutoff introduced 
by N means that the timescale is not resolved below 

r = T/N (5) 

and thus a measurement of h{t) yields effectively N data 
points. 

As one can see from Eqs. Q and ©, the amplitudes of 
the Fourier modes are independent, Gaussian distributed 
variables - but they are not identically distributed. In- 
deed, the fluctuations increase with decreasing wavenum- 
ber, with power spectrum 

(ICnP) CX i-, (6) 
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FIG. 1: Typical profiles of various 1//" signals of length N = 
T/t = 8192. Note that, contrary to the visual illusion, the 
< a < 1 surfaces are flat, while the 1 < a < oo signals are 
rough. In the former case, the amplitude of the signals is size 
independent, while in the latter case the amplitude diverges 
with system size. For ease of comparison, we have rescaled 
the signals to be approximately equal in height. 



as befitting a 1//" signal. 

By scanning through a, systems of wide interest 
may be generated. For example, a = 0, 1, 2, 4 cor- 
respond respectively to white- noise, 1//- noise |23, an 
Edwards- Wilkinson interface |29l or Brownian curve, and 
a MuUins-Herring interface pBOl l3l| or random accelera- 
tion process [s^. 

An important feature of 1//" signals is that correla- 
tions may be tuned by the parameter a. Indeed, as one 
can see in Fig^ an a-scan leads us from the absence of 
correlations (a = 0, white noise) to the limit of a deter- 
ministic signal {a = oo). In between a — and a — > cx), 
ttc = 1 separates decaying (0 < a < 1) and strongly- 
correlated (1 < a < oo) signals. As we shall see, the 
extreme statistics is different in these two regions, thus 
it may be worth spelling out the distinctions between de- 
caying and strong correlations. We therefore briefly de- 
scribe some known results regarding the correlations in 
1//" signals that will be relevant to the understanding 
of the rest of the paper. 

A global (integral) characteristic of correlations is 
given by the mean-square fluctuations of the signal, called 
roughness or width in surface terminology |33ll34| 



N/2 



W2 



[hit)-h]^ = 2j2 



(7) 



n=l 



second equality shows that W2 is the integrated power- 
spectrum of the system. This quantity has been much 
investigated [23, and its probability distribution will 
be compared to the extreme statistics of the surface in 
Sec. IXII For the present purpose it is sufficient to recall 
that the ensemble average over surfaces, (W2), yields the 
following asymptote for large system sizes (T 00) 




for I < a < 00 

for a = I 

for < a < 1 . 



(8) 



Thus the fluctuations diverge with system size for 1 < 
a < 00 in contrast to the finite fiuctuations in the < 
a < 1 regime. Since diverging fluctuations are the sign 
of strong correlations, this gives a reason for separating 
the < a < 1 and 1 < a < 00 regions and attaching 
the name of decaying and strong correlations to each, 
respectively. 

A more detailed characterization of the a-dependence 
of the correlation can be obtained by examining the cor- 
relation function Ca{t,T) = {h{t')h{t' + t)) itself. A 
simple calculation shows that the limit T 00 and 
t/T ~* finite yields the following scaling form 



Ca.{t,T)^T"-'Fa,{t/T), 



(9) 



and that the nature of the correlations follows from the 
properties of scaling function Fa- 

For 1 < a < 00, the scaling function is of order 0(1) 
and Fa{t/T ^ 0) is finite. As a consequence. 



c„>i(i,T)^r"- 



(10) 



so that the correlations diverge in the T 00 limit. 
The divergence is also present for a = 1 but it is only 
logarithmic, Ci{t,T) ~ In (T/r). Systems with 1 < a < 
00 can therefore be regarded as strongly correlated. 
For < a < 1, the correlations are 0(1) since the 



scaling function behaves as Fa{u) 



for u ^ 1 and. 



consequently, one has a power law decay of correlations, 
independent of system size 



Ca<i(i,T)^lAi" 



(11) 



where the overbar indicates an average over i, and the 



In the bulk (u ~ 1/2), the correlations quickly approach 
zero, Ca ~ i/T^-o' in the T ^ 00 limit. The correlations 
disappear entirely for a = since, in this case, h(t) are 
i.i.d. variables. Systems with < a < 1 have decaying 
correlations hence the name used for their identification. 

Thus we see how the regions < a < 1 and 1 < 
a < 00 are distinguished. Furthermore, we also have a 
characterization of correlations taken into account when 
we study the EVS of periodic Gaussian 1//" signals. 



III. EXTREME STATISTICS: TECHNICALITIES 

The quantity of interest is the distribution function 
P{hjn) of the maximum height h„i of the signal mea- 
sured from the average, as defined in Eq. In order to 
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construct the histogram for the frequency distribution of 
h^, we generate a large number (« 10^ — 10'') of signals, 
as prescribed by the action S[ck;Q:] in Eq. Q. Each 
signal is Fourier transformed and the real-space signal, 
which has zero average (cq =0), is used to determine the 
value of hm- Finally, the hm-s are binned to build the 
histogram for the MRH distribution. 

Since hm is selected as the largest from N — T/t num- 
bers, P{hm) obtained by the above recipe depends on N. 
The goal of EVS is to find the limiting distribution which 
emerges for N oo 

P(z) = lim a„P„(/i„ = a„z + 6„) . (12) 

Here ajv and biq are introduced to take care of the pos- 
sible singularities in {hm)N and in cr^ = {(hm — {hm))"^) 
(one expects e.g. that {hm) n^oo oo for distributions 
with no finite upper endpoint). 

For any finite iV, the parameters a„ and 6„ can be 
related to {hm) ^ and ct„ and, in practice, one builds a 
scaled distribution function where a„ and 6„ do not play 
any role. In the following, we shall employ two distinct 
scaling procedures. If the large N behaviors of {hm) „ 
and (T„ coincide (e.g. {hm) ^ ~ ""n ~ then we use 
scaling by the average by introducing the variable 

X = hm/{hm)„ (13) 

which ensures that {x) — 1 and makes the corresponding 
scaling function 

$(2;) = lim ^^{x) = lim {hm) ^P„{{hm) ^x) (14) 

devoid of any fitting parameters. 

If {hm) M and cr„ scale differently in the large N limit 
then the above procedure leads either to a delta func- 
tion or to an ever widening distribution. One can deal 
with this problem by measuring hm from {km) in units 
of the standard deviation i.e. by introducing the scaling 
variable 

Using y will be called a-scaling and the corresponding 
scaling function will be denoted by $(y). Provided the 
limit N ^ oo exists, 

$(y) = lim l>„(y) = lim ct„P„ ((/i™)„ + (j„y) (16) 
is again a function without any fitting parameters. 



IV. EVS IN THE REGIME OF DECAYING 
CORRELATIONS (0 < a < 1) 

In the white-noise limit a — > 0, each point on the sig- 
nal constitutes a random i.i.d. variable with Gaussian 



distribution. Under these conditions, the MRH limit- 
ing distribution falls under the domain of attraction of 
the Fisher- Tippett-Gumbel distribution P, In fact, 
in the range < a < 1, it has been shown that the 
decaying correlations are too weak to change the FTG 
limit |26| . Therefore, in the regime of decaying correla- 
tions, the MRH statistics of 1//" signals may be said to 
be universal. 




-2 2 4 6 

y 

FIG. 2: Numerically constructed MRH distributions for a — 
0.5, 0.8 and 1.0 for system size N = T/t = 16384. The FTG 
distribution is shown by the thick black line. Each distribu- 
tion is rescaled to zero mean and unit standard deviation. It 
should be mentioned that the q < 0.5 curves are not displayed 
since they are indistinguishable from the a — 0.5 case. 

However, in the case of i.i.d. random variables drawn 
from a Gaussian parent distribution (i.e. for a — 0), 
it has also been established that the convergence in N 
towards the limiting FTG distribution is logarithmically 
slow P,inEin3- Therefore, in practice, the MRH distri- 
bution may appear different from FTG. An even worse 
rate of convergence may be expected with increasing a, 
since, heuristically, increasing correlations decrease the 
effective number of degrees of freedom. In Fig. |21 we 
illustrate this trend by comparing numerical MRH dis- 
tributions for a range of a but fixed N with the FTG 
limiting distribution. For a < 0.5 the numerical distri- 
butions are practically indistinguishable from the case 
a — 0.5. This figure serves as a warning when comparing 
real-world data with known extreme value distributions. 

We note here that Eichner et al. [s^ have recently 
investigated EVS for a = 0.6. Although they do not 
spell it out explicitly, their Figs. 2 and 3 do demonstrate 
that the convergence at a = 0.6 (Fig. 3) is slower than at 
a — (Fig. 2), in agreement with our findings described 
above. 

In order to shed more light on convergence rates to- 
wards limiting distributions, we have measured the skew- 
ness 7i = K-^j k^^"^ where k„ is the n-th cumulant of the 
MRH distribution function. The results for a range of a 
and N are displayed in Fig. |31 From this plot one can 
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discern a number of remarkable features. First, in the 
range < a < 1, we note that the measured skewnesses 
are far from the skewness of the FTG distribution (ap- 
proximately 1.140...), even for the largest system size 
available. Second, for a — 0, we know theoretically that 
the convergence rate is logarithmically slow, but, some- 
what surprisingly, this convergence rate appears to be 
shared for all a < 0.5, after which convergence slows 
down markedly. Thus, the universality in the ultimate 
limiting distribution for < a < 1 may not carry over 
to a universality in the finite-size corrections. Note that 
if we did not know the limit but would try to determine 
it from finite-N skewnesses, then for a ^ 1 we would be 
wrong to conclude that the asymptotic value had nearly 
been reached. 
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FIG. 3: Strong finite-size effects at low values of a as seen in 
the scaled skewness, 71 = K-ilr^J"^ of the MRH distribution. 
Note the exponential increasing system sizes N = T/t used 
for comparisons. Note also that the scale of the horizontal 
axis changes at a = 2.5. The limiting value of 71 for a < 1 
is 71 ~ 1.140. Other exactly known values are 71(0 = 2) ~ 
0.700 and 71(0 -^00)^ 0.631. 



The case of strong correlations is discussed in the fol- 
lowing sections. Here, we just observe that the skew- 
nesses for a > 1 rapidly collapse for different iV, and 
that they are virtually indistinguishable from each other 
for a > 1.5. In this case we may be quite sure that the 
skewnesses have practically reached their limiting values, 
since they match their corresponding theoretical values 
for a = 2 and 00 with high accuracy. As we shall argue 
in Section lVIIII in contrast to the very slow convergence 
for < a < 1, convergence rate improves as it becomes 
a power law for a > 1. 



V. STRONG-CORRELATION REGIME: EXACT 
RESULTS FOR a = 2 AND a = 00. 

The 1 < a < 00 region is characterized by diverging 
mean-square fluctuations (see Eq. ijSJl). Since {hm) > 



(^2), this is also a range where the characteristic scale 
of {hm) diverges with the size of the system at least 
as {hm) T^""^)/^. An important exact result in the 
strongly-correlated regime is related to the Brownian ran- 
dom walk (a — 2). Majumdar and Comtct 7, 8] have 
shown that (hm) VTand, furthermore, they calculated 
the MRH distribution using path-integral techniques as 
well as by making a mapping to the problem of the area 
distribution under a Brownian excursion [s^ The 
resulting distribution is known as the Airy distribution. 
Under scaling by the average {x = /im/(^m))j the Airy 
distribution can be written as follows (note that slightly 
different scaling has been used in 0,|1|) 



$(2;) 



8%/3 

/^3,10/3 



(17) 

Here U (a, 6, z) is the confluent hypergeometric function 
and Vn = (2/7r) • (2a„/3)'^ is related to the n-th zero —an 
of the Airy function. 

The small and large x asymptotes of ^{x) have also 
been calculated 0, with the results 



$(a: 0) 



8\/3wi 



3/2 



and 



In $(a; — > 00 



(18) 



(19) 



It is noteworthy that the above asymptotes are quite close 
in functional form to those obtained for the width distri- 
bution of the Edwards- Wilkinson model 

The plot of ^{x) is shown on Fig. 0^ where a rather 
fast convergence to the limiting function can be seen (the 
convergence rate is T~^/^ as calculated in ^3). It is re- 
markable that the convergence is even faster if a-scaling is 
used (Fig.^lD). The reason for this is the finite-size scal- 
ing of higher cumulants of the MRH distribution func- 
tion. At this point we present this just as a numerical 
observation. A detailed study of the finite-size scaling of 
MRH will be pubHshed separately [i^ . 

The review of the properties of the MRH distribution 
for a = 2 presented above gives a guidance for the discus- 
sion of EVS in the strong correlation regime. As we shall 
see below, the basic properties of EVS {{km) ^ \J (^2), 
the general shape of $(a;), the structure of the small- and 
large- a; asymptotics of $(a;), the fast convergence to the 
limiting function) are similar in the whole 1 < a < 00 
region. 

The other analytically solvable case is the a 00 
limit. Indeed, here only the n = \ mode survives, 
and the resulting signal h(t) = |ci | sin (27rt/T -f 93). 
Consequently, h^n = \ci\ and the distribution of /i,„ 
is just the distribution of |ci| given by P(|ci|) ~ 
|ci| exp [— (27r)°r^^"|ci p]. Using the average scaling^ the 



scaling function <I>^ 



{hm)P{{hm)x) becomes 



-x9{x)e 



-TTX /4 



(20) 
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FIG. 4: Convergence of the MRH distribution to its exactly 
known N ^ oo limit 0. Results for same system sizes are 
displayed in both panels using scaling by the average {x = 
hm/ihm)) and cr-scaling [y = {hm — {hm))/a\ in the upper 
and lower panels, respectively. Very small system sizes were 
also displayed in order to demonstrate the remarkably fast 
convergence by cr-scaling. 



were 9{x) is Heaviside's step function. Comparing the 
above expression with the asymptotes (|18|) and (|19|l , one 
can see that, in addition to the disappearance of the small 
X singularity, the large x asymptote has also changed by 
an extra x factor. 



MRH distribution function (Sec. IVII|I from which rather 
general conclusions can be drawn and, furthermore, as 
indicated by the simulations, the results can be extended 
to any 1 < a < oo. 

The construction of the MRH distribution function 
in the path integral approach involves the calculation 
of a normalization factor which, in turn, requires the 
knowledge of the propagator (also called two-time Green- 
function, or transition probability) of the random acceler- 
ation process. Here we compute this propagator, i.e. the 
probability density of a position of the stochastic path at 
some time t conditioned on the initial point. 

The equation of motion of the 1//^" trajectory reads 



{t) 



d"hit) 

di" 



(21) 



where ^(t) is white noise with zero mean and correlation 
= S{t - t'). Note that h{t) corresponds to the 
n — 1-st integral of the Brownian random walk trajectory. 
Eq. I|21|) can be rewritten as a vector Langevin equation 



Zk = Zk-1- 



(fc = 2,3,...n), zn = h. (22) 



For n = 1 we have the usual random walk, for n = 2 the 
random acceleration problem [s^ , also extensively stud- 
ied, while for higher n-s one can speak about the gen- 
eralized random acceleration processes |43l . We are 
interested in the conditional probability that after time t 
the trajectory is at 2; = (zi, 2:2, ... , Zn) provided it started 
from . In the following, we denote this propagator by 
Gn{z\z'^;t). Its subscript indicates the dimension of the 
vector arguments, and it obviously satisfies the recursion 
relation 



Gn{z\z^;t)^ J dz„+iG„+i(2N°;i) 



(23) 



where the integration eliminates the dependence on 
too. The propagator has Dirac delta initial condition, 
Gn{z\z°;0) = (5(")(2 - 2°), and satisfies the Fokker- 
Planck equation, obtained in a standard way from the 
Langevin equation , 



dfGn — —H^Gn 
1 



n-l 
k=l 



Zkdi 



k+l, 



(24) 
(25) 



VI. THE PROPAGATOR OF THE 
GENERALIZED RANDOM ACCELERATION 
PROCESS: a = 2n 

The derivation of subsequent analytical results on the 
scale of (hm), and on the small hm asymptote of the 
MRH distribution is based on the observation that 1//" 
signals are actually paths of generalized random accel- 
eration processes, provided that a = 2n is an even in- 
teger. This allows a path-integral representation of the 



where dt and dk are derivatives with respect to t and 
Zk, respectively. The superscript of refers to the fact 
that we consider here the time evolution (|21() without 
further constraints. The propagator has been calculated 
in previous studies up to n = 5 0, . We make an 
ansatz that matches these functions and we show it to 
be valid for general n 



Gn{z\z";t) = II g{a' 



;o-fc), 



(26) 



fe=i 
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where Q{z;a) = exp (— z^/2cr^)/ y/2na is a Gaussian PDF 
with zero mean and variance tr^, z'^ is the initial condition 
vector, and the vector a'' only has nonzero components 
for i = 1,2, . . . ,k, i.e., for I > k we have af — 0. In order 
to remove ambiguity we set — 1. The a*^, a^'^"', (jfe are 
time dependent quantities to be determined. The above 
formula amounts to the recursion relation 

Gn{z\z";t) - Gn-i{z\z";t) 

xg(a"-;z-aO'"-;zO;f7„). (27) 

Substitution of this ansatz into 124|) leads to equations for 
the unknown parameters. The solution of the equations 
as described in Appendix ^ yields 



(-t)"-'=(n + A:-2)! 
2"-i(2n-3)!!(A:-l)!(7i-A:)!' 



-fe. 



i(2n-l)/2 



2"-V2n- l(2n- 3)!! 



(28a) 
(28b) 
(28c) 



For illustration, we use the above expressions to calcu- 
late and display explicitly the n = 4 propagator. Noting 
that the original coordinate, velocity, acceleration, and 
its time derivative are given by 

h — Z4, V — Z3, a — Z2, a = zi, (29) 

respectively, the propagator can be written in the form 

720\/T05 



G4(z|z°;t) 



-A/2 



(30) 



where — A/2 is the sum of the exponents of the Gaussians 
in H2t)|) . namely 



k=l 



with 

A2 
A^ 
Aa 
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^3-Z§--(z2 + Z2") + — (Zi-Z?) 



(31) 



(32a) 

(32b) 
2 

(32c) 



Zi-Zl- -(Z3 + Z3") + — (Z2 - Zl) , 



(32d) 



Up to the k = Z term this incorporates the propagators 
of the random walk, fc = 1, random acceleration, k < 2, 
and random velocity of acceleration, fc < 3, and the above 
expressions are in agreement with previous results [44j . 



Note that, independently of k we have a* 
and 



0,k 

H-1 



t/2, 



(33) 



but for / < fc — 2 the af-s will vary with both I and fc. 

Later, for the construction of the formula for the MRH 
distribution, we will need a special property of the prop- 
agator. Namely, if we consider the propagator of a peri- 
odic path of length T and integrate it over the common 
values of the velocity, acceleration, etc., at the endpoints, 
we get the surprisingly simple result 

/n— 1 
ndzOG„(z°|;z°;T)-T-("-i/2)(27r)-i/^ (34) 
fe=i 

Indeed, the periodic propagator does not depend on 
the integration over z^_-y cancels the normalizing con- 
stant of the n-th Gaussian but brings in a factor of 1/T. 
The integration over zJJ„2 does the same with the n — 1- 
st Gaussian, and so on, until finally we are left with the 
norm factor of the n = 1 Gaussian, 1/^/2ttT, divided by 
T""^, as shown in H34|) . The key to this remarkable can- 
cellation of the total numeric prefactor of the propagator 
is that H33|) holds uniformly for all fc-s. 



VII. PATH INTEGRAL FORMALISM AND THE 
TRACE FORMULA FOR THE MRH (a = 2n) 

For a = 2 Majumdar and Comtet 0, ll| introduced 
a path integral representation of the MRH distribution. 
The technique allowed for the formulation of the PDF 
in terms of the spectrum of a quantum mechanical, one- 
dimensional, Hamiltonian H with a hard wall and else- 
where linear potential, through the trace of e^^^, valid 
in the case of periodic boundary conditions. The spec- 
trum is known to consist of the Airy zeros, so the trace 
formula resulted in the PDF called the Airy distribution. 

In what follows we show that, in the case of periodic 
boundary conditions, for a general a = 2n, n = 2,3, . . . 
an analogous trace formula holds. Remarkably, the for- 
mula turns out to be essentially the same as in the a = 2 
case, with the only difference that now a generalized 
"Hamiltonian" Jf„ appears. However, the ff„, a differ- 
ential operator in an n dimensional space, is no longer 
Hermitian. Whereas we shall not solve the spectral prob- 
lem necessary for the calculation of the MRH distribu- 
tion, this formulation will allow us to (i) determine the 
scale of the MRH as function of T, and (ii) give explicitly 
the initial asymptote of the PDF, with the only undeter- 
mined parameter being the ground state energy of the 
Hamiltonian ff„. What is more, the results (i-ii) will 
lend themselves to a continuation to real as, so the use 
of the path integral technique extends beyond its origi- 
nal region of validity, the generalized random acceleration 
problem a = 2n. 
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We begin with the probabihty functional of a periodic 
path h(t), where h is measured from the time average, 



P[h{t)] = Aexp(^~^J^ dt[h\ 



n] ^,m2 



it)Y 



xS / dth{t) 



(35) 



Following 0, Q we have introduced a normalizing coeffi- 
cient A, ensuring 



where the mark PBC:2;*^ refers to the time derivatives 
at the ends fixed at /i['=l(0) = /iW(r) = z°_^, k = 
0, . . . , 71 — 1. The PpBci^'^) is the joint probability den- 
sity of h'^'^^t) = z^_,^-s in a periodic path at any fixed 
time t, so as a byproduct we obtained that density in 
terms of the propagator, explicitly given in Sec. I VII That 
joint probability density is obviously normalized to unity. 
However, we know from Eq. (|34|l that the integral of the 
r.h.s. is independent of the z^^j^-st variable, and there- 
fore we have 



A 



'2tt. 



(41) 



V^h{t)P[h{t)] 

PBC 



(36) 



where PBC indicates that periodic boundary conditions 
for all derivatives of the path up to is understood. 

The measure 2?„ft.(t) is defined such that the propagator 
Gn{z\z'^]T) of Sec. I VII is a path integral without extra 
normalization, and the boundary conditions of the inte- 
gral are specified by the arguments of G„, i.e. 

G„(^|^"; T) = jn^h{t) exp £ dt [/jN (t)]^^ ,(37) 

with /i["-'^l(0) zl, /i["-'=l(T) = Zk, where k ^ l,...n. 
This 'Dnh{t) is in fact the measure leading naturally to 
the quantum-mechanical-like operator representation of 
the path integral 



For n — 1 the normalizing coefficient derived in 0, Q is 
recovered. 

The integrated distribution M{h^\ T) of the MRH, i.e., 
the probability that the maximum does not exceed /im, 
has been formulated in terms of a path integral in [3,l3l- 
That expression is valid for any path density P[h{t)\ and 
reads formally as 



M{h^-T) = 



dhP{h;T) 

G 

Vh{t)P[h{t)]l[e{h^-h{t)).{42) 



PBC 



G„(2|z°;r) 



z\e-"^^\z° 



(38) 



where is given by Eq. (HSJ), and \z") {{z\) are its 
right (left) eigenvectors corresponding to the n dimen- 
sional positions indicated therein. Note that since this 
Hamiltonian is non-Hermitian for n > 2, the left and 
right eigenfunctions are different in general. 

A third version of the propagator we shall utilize comes 
from a path integral by a measure of one order lower as 



Note that here P(/im; T) is the density of MRH. Changing 
the integration variable and then introducing the hard 
wall potential Vo{h) — oo for h < Q and Vo{h) = for 
h > 0, one obtains 

M{K,-T) = [vh{t)P[K,-h{t)]\\e{h{t)) 

JPBC ( 

^Vh{t)P[K, ~ /i(t)]e-^o'''^*^'(''(*».(43) 

PBC 

Using the specific form (|35|l of the probability functional 
we find 



G„+i(z|z";T) = JVnh{t)exp^~^J^ dt [/i^ (t)]2^ M(/i„,;T) - A J^Vh{t) S (^K,T - dt hit)^ 



x5|^z„+i-y /i(t)dt - j .(39) 

Here the Dirac delta produces the normalized density for 
the added area variable Zn + 1- Note that if we take 
into account Eq. (|37|l then the consistency relation 123() 
immediately follows. 

In order to determine the normalization coefficient A 
in (jSni, we express the equal-points propagator comple- 
mented with the area variable set to zero at both ends. 
By integrating over the path except for a single point and 
using Eqs. and 

Pfbc{z°) = I VMt)P[Ht)] 

= AGn+iiz°,0\z°,0;T). (40) 



X exp 



T 
dt 



i(/iN)2 + y„(/,(t)) 



44) 



Next, we introduce the scaled Laplace transform of the 
integrated MRH distribution 



K{u;T) = T f 
Jo 



dh„ 



-uh„,T 



M(h^-T) 



A I Vh{t) 

'PBC 



X exp 



dt {'^{h^'^^f +uV{h{t)) 



(45) 



where the potential V{h) = oo for /i < and V{h) = h 
for /i > 0. 
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In order to find K{u; T), we write down the evolution 
equation for the PDF of the position and its derivatives 
Pn{z;t) corresponding to the above path probability 



(46) 
(47) 



where -ff^ was given in (|25() and the variables Zk are de- 
fined by 122|l . Thus the Laplace transform can be written 
in short as 



(48) 



K{u;T) = ATiexp [-Hniu)T 



It is straightforward to show that the eigenvalues 
En,Lj{u) of Hn{u), where lo summarizes all discrete in- 
dices, obey a simple scaling in u. For that purpose, let 
us consider the eigenvalue problem for /i = z„ > 



-dl + Zkdk+i 



k=l 



^ = £;„(u)^(49) 



and apply a scale transformation by substituting 

Zk, u^En E^. (50) 



u'^'^Zk 



We recover an equation free of u, if all powers multiplying 
various terms are the same, that is 

- 2/3i = /3i - /32 = • • • = Pn-i - /3„ = /3„ + 1 = <5, (51) 

so f3k ^ {n — k + 1)5 — 1 and also — 2/3i = S. Hence S = 

2;^, so the eigenvalues scale like En,u{u) — en.ujU^^ , 

where en,uj is the spectrum of Hn = Hn{l)- 

It thus follows that, using (|41|l . we get the scaling re- 
lation for the Laplace transform of the integrated distri- 
bution 



K{u;T) = r"+2/C(ur"+2) 

/C(s) = \/2^Trcxp (^-7?„s^ 

= V2Tr exp {-en^LjS^ 



(52) 



(53) 



Hence, using (gSJ, we obtain for the PDF of the MRH 
P{hi^; T) and its moment generating function G{u] T) in 
scaling forms 

P{h^;T) = 9^„M(/i,„;r)=T5-"P(/i,„r^-»), (54) 

G{v-T) = ( dh,^P{hn,;T)e-"''- ^G{vr'-^),{55) 
Jo 



where 



G(s) 



oc 



dz P{z)e 



'0 



iS2„ + l 



27r s Tr exp (^—Hn 
V2tt s ^ exp (^-en,ujS^ 



(56) 



Note that the same symbols P, G are used for single- and 
double-argument functions, but that should not cause 
confusion. Remarkably, the trace formula is exactly the 
same as in the case of the simple random walk n = 1, 
with the Hamiltonian Hi replaced by iJ„. Note that, 
in the special case of n = 1, the scaling function of the 
MRH distribution l|17l) is ultimately recovered from the 
above trace formula 8]. 

The scaled moment generating function Eq. to- 
gether with the preceding scaling formulas are our main 
result here. In the next two sections, we shall exploit the 
above results to draw conclusions about the scale and 
the small argument asymptote of the MRH distribution 
function. 



VIII. STRONG-CORRELATION REGIME 

(1 < a < 00) 

In order to evaluate the trace formula one would need 
the energy eigenvalues of iJ„(l). Although they are 
known .7] only for n — I, assuming that these eigen- 
values exist, the scale of the MRH in T can be derived 
since Eq. (|54|) yields 



(57) 



As one can see, the scale of (hm) is the same as that 
of the square root of the roughness i.e. we have 
(/^m) ^ \J {102) just as in the case of random walks 
It should be emphasized that while the above reasoning 
holds strictly for a = 2n, the exponent can, in fact, be 
continued naturally to real values. Thus it is plausible to 
surmise that is the scale of the MRH for any a > 1. 

This power emerges quite sharply for a — 1.3,1.6,2,4 
in numerical simulations as shown on the first panel of 
Fig. El 

Fig. El also displays the 2nd and 3rd cumulants of 
P{hrn,T). We can observe the emergence of well defined 
scaling with T 



(58) 



The scaling exponents are again equal to those of the 
cumulants of the width-distribution provided the hm ~ 
y/w2 correspondence is used. This suggest that there is 
an intimate connection between the fluctuations of MRH 
and those of the signal width. 

In order to see how the general shape changes as a is 
increased, we have performed simulations as described in 
Sec, mil The results are shown in Fig. El where we used 
scaling by the average to present the scaling functions 

The main features can be readily seen. The scaling 
function is a unimodal (single peaked) function which 
spreads out as a increases and approaches its a — *■ 00 
limit (see Eq. (|20|l ') rather fast. This is not entirely sur- 
prising since a glance at Fig. ^ convinces one that the 
a = 10 signal already consists of a single mode for all 
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N 

FIG. 5: Cumulants, Kk, of the MRH distribution for various 
a > 1 showing that Kfe(a) scales with system size, A'^ — T/t, 
as Kfe(a) ~ TV*^^""^-*/^. The straight lines have the appropriate 
asymptotic slopes k{a — l)/2. 
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FIG. 6; Numerically constructed MRH distributions for vari- 
ous 1 < a < oo. The a = oo curve is the analytic result given 
m Eq. ll^ni. System sizes A'^ < 16384 were sufficiently large 
to observe the convergence of the PDF-s within the width 
of the lines drawn. Note that the a — S results are almost 
indistinguishable form the a oo limit. 



The large x limit is harder to treat analytically and we 
have only numerical evidence (Fig.CJ) that the asymptotic 
behavior for large x is given by 

$(x) ~ Ca;^e-^="' (59) 

where the parameters B, C, and 7 depend on a. The 
above functional form is consistent with the exact re- 
sult at a = 00 (see Eq. (HOI)- At a = 2, the ansatz 
of a Gaussian decay was shown to be in agreement 
with the large-order moments of the distribution func- 
tion. However, the possibility of a prefactor x'^ was not 
excluded by the analysis. We found that the generalized 
asymptote with 7 « 2 gives a superior fit to the 
large- a; {x > 1.5) behavior of the exactly known PDF. 

We have also fitted our numerical data in the region 
X > 1.5 for larger as, resulting in 7 « 1.4 and 7 w 1.1 
for a = 3 and 4, respectively. The general trend of the 
exponent 7 with increasing a is consistent with the a — > 
00 limit of 7oo ~ 1- 



IX. INITIAL ASYMPTOTE 



practical purposes, and thus the MRH distribution will 
be very well approximated by the $00 (2;) function. 

The function decays to zero extremly fast in the x — > 
limit. The nonanalytical behavior and the actual func- 
tional form at small x will be the subject of the next 
Section. Here, we call the reader's attention to the fact 
that the region where the nonanalytic asymptotic behav- 
ior dominates is shrinking as a increases and, according 
to Eq. H20|) . entirely disappears in the a 00 limit. 



The trace formula l|5t)|) allows us to perform an asymp- 
totic analysis of the MRH distribution for small argu- 
ments. The calculation is based on the large s behavior 
of the moment generating function (|56|l , wherein we as- 
sume that there is a positive, a — 2n-dependent, non- 
degenerate ground state energy 60(0;), which gives the 
leading term of the sum (while a is strictly even, several 
results will lend themselves to continuation). Under this 

assumption, the PDF in the scaled variable z = hraT^~ 
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is asymptotically given by 
ds 



P{z) 



ds r— 

sv2tt exp sz 



(60) 



The above integral can be calculated using the saddle 
point method. For small z, the saddle point of the expo- 
nent is located on the real axis at 



2eo 



l)z 



(61) 



and the integral in the neighborhood of the saddle point 
reduces to evaluating a Gaussian integral which yields 
the following asymptote 



P{z) w Cz-^exp(-B/z'^) , 
where the parameters are given by 
2 



B 



C 



a - 
2a 



1 



1 



a 

a - 



2eo 



a-1 



1-1 
2eo 



1 



(62) 

(63a) 
(63b) 

(63c) 
(63d) 



One should note that the exponents (3 and 7 depend only 
on a while the amplitudes also depend on the ground 
state energy, eo(a). The value of eo{a) is known only for 
a = 2 where eo(2) = ai/v^, with ai — 2.3381 being the 
absolute value of the first zero of the Airy function [3, |^ . 

It should be emphasized that wc did not scale the mean 
to 1, being ignorant about the full PDF as well as its 
mean for general a. So if comparing the above formula 
to the numerically scaled PDF as function oi x = z/ {z) 
then the factors B, C will change and become fitting pa- 
rameters. Figure demonstrates the fit of (|62|l to sim- 
ulation results for several a-s, and we find that the fits 
are excellent in a surprisingly large interval. It should 
be noted that in the large a limit the initial slope is 
positive, so one expects a decreasing range of validity of 
the asymptote for increasing a, nevertheless, the fit on 
Fig.[7|is quite good even for the largest a. The case a = 3 
demonstrates the continuation of the a = 2n based for- 
mula, and suggests that naive continuation of at least the 
exponents /3, 7 in l|63|) is justified. 

Returning to the problem of scale-dependence of the 
amplitudes B and C, we note that even if the full PDF is 
unknown, one can construct a parameter from the small- 
X asymptote which does not depend on the scale. In order 
to see this, let us consider scaling by the average. With 
the rescaled variable x, one has the PDF as (z) P{x (z)) 




FIG. 7: MRH distributions for (a): a = 2, (b): a = 3, 
and (c): a = 4, calculated for system size N = 16384 (solid 
lines). The small and large x asymptotes (dashed lines) are 
also shown. The small x behavior in the range < a; < 1 
is fitted to the functional form 1621 1 with the exponents /3, 7 
taken from I63II . The prefactors B,C are fitting parameters 
(note that the formulas in (I63II contain an unkown parameter 
eo). Large x data in the range a; > 1.5 are fitted to the form 
15911 ■ where B, C and 7 are fitting parameters. 



and writing it again in the form (|62|l yields the following 
change of the amplitudes 



B' 



B 



c 



c 



\7-l 



(64) 
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It follows from the above expressions that the following 
combination 

^ ^ C ^ C 

remains independent of any scale change. 

We should reiterate that the energy parameter eo(a) 
is not known generally, but it is plausible to assume that 
it is a well defined number. It may be determined nu- 
merically for a = 2?T, by a direct study of the correspond- 
ing local Hamiltonian. Remarkably, however, the above 
asymptotic formula allows for the computation of eo{a) 
for any a > 1 from a numerical fit of the simulation 
result. Thus, precise MRH statistics effectively extract 
the ground state energy level of the Hamiltonian without 
solving the corresponding differential equation. Contin- 
uation of H62|l for a ^ 2n is also natural here, but in 
this case we have a non-local Hamiltonian, whose spec- 
tral problem would be an even more challenging task to 
solve. Unfortunately, very high precision simulations are 
required to determine the ground state energy from the 
small- a; asymptote. In particular, our simulated data did 
not even allow the computation of the ground state en- 
ergy to within a factor of 2 for the case of a = 2 where 
the lowest eigenvalue is known. 

X. MRH DISTRIBUTION FOR LARGE a 

We have calculated the MRH distribution for the a 
oo limit in Sec. O There we found that only the n = 1 
mode survives and, as a result, the PDF H2U|) emerges. 
Here we discuss a procedure for perturbatively computing 
the leading corrections to H2U|) by taking into account the 
modes n = 2, 3, . . . . 

First we reiterate that the amplitude of modes c„ obey 
the distribution with action ^ and measure proportional 
t'-' Yin ^(cn)cndc„. Thus, Separating the n = 1 mode, the 
path in Fourier representation is written as 

oo 

h{t) — ai sm(t) + SnCLn sin(nt + ipn), (66) 

n=2 

where the e„ = 1/n"/^ is the mean square root deviation 
of the amplitude of the n-th mode, and the a„ — Cn/SnS 
are i.i.d. variables distributed according to 

Po{z) ^ 2z0{z) exp(-z2) . (67) 

Finally the phases (/?„ are independent and uniformly dis- 
tributed in [0, 27r]. The n = phase is omitted, because 
the choice of the origin is arbitrary. Obviously, x mea- 
sures the height from the time average of the path, which 
is here set to zero. Note that now time t is in units of 
2tt/T. 



The leading correction from higher frequency modes 
can be calculated independently for each mode, thus here 
we only consider the n-th mode. Then the path is 

h{t) ~ ai sin(t) + e„a„ sm{nt + ipn) (68) 

and the calculation to leading order is straightforward. 
We compute the maximum of the path and then, know- 
ing the distribution of all parameters therein, we can de- 
termine the PDF of the maximum. The details are pre- 
sented in Appendix ^ where we obtain the perturbed 
PDF for h,nT^ = z as 

P{z) ^ Po{z) + elP^Az) (69) 

with 

= {l-ny2)5{z) + z6'{z)/2 

+e-^'6'(z)(2z3 + (ri2-3)z). (70) 

The singular part needs some explanation here. As has 
been discussed in Sec. IIXI for finite a-s the PDF starts 
nonanalytically with zero initial slope for finite a-s, in 
contrast to the a — 00 case, where the PDF has a fi- 
nite slope. The nonanalyticity is not expected to be re- 
covered by any expansion. Nonetheless, the formal ex- 
pansion gives an explicit correction function P2,nj with 
delta-singularity at the origin. It is plausible to conclude 
that while the expansion cannot be correct overall, the 
singularity "tries" to take care of the nonanalytic differ- 
ence in the small-z behavior, while the nonsingular part 
is expected to be a faithful correction for z > 0. This 
leaves open the possibility that the large a expansion is 
not convergent, rather it is asymptotic. 

Next we scale the PDF to unit average. Using the 
result l|B8l) from Appendix IbI one finds 

<^{x)^<^>^ix)+el<P2A^) (71) 
where ^^{x) is given in H2U|I and 

*2,„(x) = {l-n^/2)6{x) + xS'{x)/2 

+e-^^'/^9{x)l{n^ - 1) (6x - ttx^) .(72) 
8 

Formally, we can sum up the leading corrections for all 
n-s. However, this is not a consistent approximation, 
because, for instance, £2 — so the leading correction 
from n = 4 is of the same order as the quadratic one from 
n = 2. Therefore we use the sum of the corrections for 
only n = 2, 3 to test the prediction. On Fig.|Hlwe display 
the correction 

A$(a;) = $(x) - $^(a;) (73) 

of the PDF from simulation for a series of a-s, together 
with the theoretical prediction for corrections added up 
from the modes n = 2, 3, 

A$(x) «e2$2,2(x)-)-£^$2,3(a^). (74) 
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FIG. 8: Difference AP(a;) between the MRH distribution for 
large a and a = oo, as defined in 1731 . We also display the 
smooth part of the leading correction 17411 from perturbation 
theory added up from the modes n — 2, 3. 



It should be mentioned that the sum of leading cor- 
rections for all modes has a prefactor X]^2('^^ ~ ^)^n — 
({a — 2) — C(Q!). This diverges for a 3+, so the series 
does not converge for a < ac = 3, which is the bor- 
derline for the differentiability of the path. We cannot 
exclude that higher order corrections, involving higher or- 
der differentiation of the path, further tighten the range 
of convergence. 



XI. COMPARISON WITH THE ROUGHNESS 
AND THE MAXIMAL INTENSITY 

Here we compare the statistical properties of hm to 
those of the roughness W2, i.e. of the mean square devi- 
ation, or width of the trajectory {T)). The latter was one 
of the first global quantities of stochastic signals whose 
scaling properties and statistics were extensively studied 
for 1//" processes pvj . 

One of the reasons for comparison is that both (hm) 
and \/(w2j scale similarly with T for a > 1 and, further- 
more, there are many common features at the level of 
their PDF-s. Namely, for algebraically diverging correla- 
tions, a > 1, after scaling by the mean, the PDF-s are 
nondegenerate (each cumulant is finite), that is, scaling 
by the mean is a natural representation of both PDF- 
s. As a ^ 1, the PDF-s scaled by the mean approach 
a Dirac-delta and, in the range a < 1, they both lend 
themselves to scaling by the standard deviation. Here an 
important difference emerges. In the range 0.5 < a < 1 
the roughness has a nontrivial PDF while below the crit- 
ical a ~ 0.5 it becomes trivial, i.e. the roughness be- 
comes Gaussian distributed. On the other hand, the 
MRH has the trivial FTG limiting distribution in the 
entire < a < 1 region. We can only speculate that the 
threshold near a — 0.5 manifests itself in the MRH distri- 
bution in its approach to the FTG limit, as suggested by 
the finite-size dependence of the simulation results shown 
in Fig.O This, however, is just a numerical observation 
without theoretical foundations as yet. 

Further motivation for a closer comparison comes from 
the similarities in the shape of the two families of PDF- 
s in the a > 1 region. First, in the a —^ oo limit, the 
PDF-s are the same if the h„i ^ ^/w2 correspondence is 
made. Second, for finite as, the unimodal PDF-s have 
asymptotes which are similar for both small and large ar- 
guments. Specifically, there is a Gaussian decay at large 
X, while the small x behavior is dominated by an expo- 
nential nonanalytic term with a power prefactor. Here, 
the comparisons can be made quantitatively for small x, 
since analytic results are available for general a. 

Last but not least, a reason for a closer comparison 
comes from the fact that the roughness can also be con- 
ceived as obeying an EVS. Bertin and Clusel 0, 
made the remarkable observation that since the rough- 
ness is essentially the integrated power spectrum, i.e., 
the sum of nonnegative Fourier intensities, it is in effect 
the maximum of positive partial sums. In general, the 
partial sums are correlated but, for the special case of 
a = 1, they correspond to the ordered sequence of i.i.d. 
variables. As a consequence, FTG distribution emerges 
for W2 at a = 1, thus providing insight to an earlier 
rather puzzling result in connection with 1// noise. 
It then becomes a rather interesting open question how 
the MRH distribution differs from the roughness distri- 
bution for a > 1 where the latter also describes the EVS 
of correlated variables. 

The initial asymptote of the MRH distribution 
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(|62|I . H63() should be compared with that for the rough- 
ness distribution obtained in appendix E of Il0| 

<i>w{x) - C^x'^-exp{~B^/x''^-), (75) 

where the parameters are given by 



Bu 

Cu 



1 



a- 1' 
3a — 1 



(a-1) 



^Q!sin(7r/a) 
(27r)^ / TT 



\/a — 1 Vsin(7r/a) 



(76a) 
(76b) 

" C(a)-^ (76c) 
(C(a)a)"5T^76d) 



with C(q^) denoting the Riemann's zeta function. Note 
that this asymptote does not contain unknown parame- 
ters such as eo in the MRH distribution. 

Interestingly, comparison with the exponents in the 
asymptote of the MRH, (|62|l . shows that the respective 
7-s are the same, if y^u^ is considered, i.e., = 7. 
Nevertheless, the respective exponents in the prefactor, 
2/3tu -|- 1 and /3, agree only at an accidental a and are 
otherwise different. 

The present results on the small- a: asymptote may also 
be compared to the asymptote of the distribution of the 
maximal Fourier intensity. It is defined as the maximal of 
the |c„p intensity components for a given realization of 
the path, which obeys some PDF if the ensemble of 1//" 
paths is considered. This was to our knowledge the first 
quantity whose EVS was studied in the context of 1//" 
signals J^J. Again, the overall shape of the PDF of the 
extremal intensity is similar to those of the MRH and the 
roughness: its initial part is suppressed nonanalytically 
and has a single maximum, before smoothly decaying for 
large arguments. There the critical ac where the FTG 
limit distribution emerges is q^c = 0, in contrast to the 
MRH and the roughness, where this critical values are 
ac = 1 and ac = 0.5, respectively. As we have shown 
in |lClj| . written with ac, the powers in the asymptotic 
formula for the maximal intensity and for the roughness 
are the same 



1 



Iw =71 = 



3{a-ac) + 2 
2(a-ac) ' 



(77) 



where the exponents /3i , 71 are defined in the same way for 
the initial asymptote of the PDF of the maximal intensity 
as Pw,7w were for the PDF of the roughness. 

In conclusion, the respective PDF-s of the MRH, the 
maximal intensity, and the roughness are similarly look- 
ing, unimodal functions, with nonanalytically slow initial 
behavior. Despite the qualitative similarities, however, 
it is clear that the three PDF-s are quantitatively differ- 
ent. This is natural since they describe different physical 
quantities. One may, however, speculate that the simi- 
lar features have their roots in the divergent correlations 
present in the a > 1 region. 



XII. FINAL REMARKS 

It should be emphasized that we are only at the first 
stages of understanding the effects of correlations on 
EVS. One of the important tasks for future studies should 
be the understanding of the convergence properties in 
the < a < 1 range. Although the limit distribution 
is known here, the convergence is extremely slow. Since 
most of the environmental time series of general interest 
(data on temperature, precipitation, etc.) correspond to 
this range, as they exhibit generically correlations with 
power-type decay, and the length of the series is naturally 
restricted, the development of a theory of finite-size cor- 
rections is important. The much discussed a — > 1 case is 
even more challenging since it appears to be outside the 
reach of present computing abilities. Thus new analyti- 
cal approaches and ideas for numerical recipes are called 
for. 

Another relevant problem is the question of bound- 
ary conditions. It is known from the a = 2 case, where 
both periodic and free boundary conditions were investi- 
gated , that the MRH distribution depends on bound- 
ary conditions. Since the analysis of a real time series 
usually means cutting it up into smaller pieces and mak- 
ing statistics out of the properties of these subsequences, 
the appropriate boundary conditions in this case are the 
so-called window boundary conditions, when the win- 
dow under consideration is embedded in a longer sig- 
nal. These boundary conditions have been discussed in 
connection with the roughness distribution of 1//" sig- 
nals [23 ■ It has been found that the limit distributions 
depends on the window size (even in the limit of large 
external system) and furthermore, the effects become 
stronger as a increases. Clearly, similar studies should 
be carried out for the EVS problem. 

Finally, it remains to be seen if the investigations of 
the effects of correlations, in particular the effects of 
strong correlations, will allow us a universal classifica- 
tion of EVS similar to that existing for thermodynamic 
critical points. 
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APPENDIX A: DERIVATION OF THE 
PROPAGATOR FOR a = 2n 

Here we show that the ansatz (|27|) indeed satisfies the 
Fokker-Planck equation H25|l with the coefficients H28|l . 
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Let us start out from (|25|l 



ri-l 



and substitute 



Gn ~ Gn-lG, 



(Al) 



(A2) 



where the arguments are understood as in H27|l . Using 
the fact that G„_i also satisfies (IA1|) . we arrive at 



Thus the full first hues on the r.h.s. of (|A7|) and HA8|I are 
equal. In the rest we change the summation variable k to 
I and then equate the respective factors of Zk and those 
of z° to obtain differential equations for the coefficients 



n-l 



.1 „i 



'1 

l=k 



• 0,n 

a I, = 



.1 



n 



'1 Zl^ 

l=k 



(AlOa) 
(AlOb) 



I n— 1 

= -dfg + diG di lnG„_i - J2 Zkdk+iG- (A3) 



fe=i 



From (EZIl we have 



dig = {alYQ" 
where, denoting a^-z — a?''^-z^ by 



(A4a) 
(A4b) 



(A5a) 



g"{x;a) = -^i^-^g'{x;a), (A5b) 



(7^ 



and, furthermore, using the full exponent of the ansatz 
we have 

n— 1 k 

9ilnG„_i =-^%(a'^--^-a"'''^2"). (A6) 

k=\ 

Differentiation of the Gaussian Q by time gives 



(A7) 



where we have condensed some z dependence by factoring 
out Q' . On the other hand, according to Eqs. l|A4p . I|A5|I . 
the r.h.s. of (|X3l) yields 



1 



dtQ = ^^{a^r ^ + 



2 ' \(Ti, 

ii y z — a -z ) 



n-l 



(A8) 



k=l 



Equating ljA7|l with IjASp should give the sought after 
equations for a, a*', a. Comparing the ^-independent fac- 
tors of g in ((X7|) and in (|A8p gives 



(A9) 



where we have used the condition that a^^ = a°'' = for 
k > I. Next, we determine the time dependence of the as 
by assuming it to be power law and requiring that terms 
in each differential equation have the same power. Thus 
we separate the time dependence, and for later purposes 
also factorize the constants as 



fn-kiii u 



(All) 



for all nonnegative integers n, k. We also set 6^ = for 
n < fc to ensure that vanishes for such indices. The 
c„-s are made unambiguous by requiring 6" 1 and, 
furthermore, since a" = 1 thus 6" = l/c„. Then (|A9|) 
gives 



'^1 = iCnf 



2n-l 



(A12) 



The parameterization in I|A11|) is justified by the fact 
that substituting it into (jAlOp the c-s disappear, so what 
remains are equations for the &-s as 



n-l 



K+i - -{n~k)hl~Y.^2l~l)h\ (Al3a) 

l^k 

Y.^2l~l)bf. (AlSb) 



°k = ^ 

" n-k 



l=k 



For a few small integer indices these equations can be 
solved, whence the following general formulas can be sur- 
mised 



hi 



hT = (-1) 



(-1) 



(n + fc-2)! 
(n- fc)!(fc- 1)!' 



n — kj-,n 



(A14a) 
(A14b) 



Note that Eqs. (EUt) and (jSISb) are homogeneous 
linear equations leaving room for overall factors in the 
solution. They are set by the conditions (i) 6" = 1 and 
(ii) = 6". Condition (i) was stated earlier below 
Eq. (|A11(I , while (ii) is equivalent to the requirement that 
Gfi depends on z„ and only through their difference. 
One can confirm proposition ljA14|l by substituting it 



16 



into (|A13|) and then using the identities 

'^i2l-l)ik + l-2y. n-l{n + k-2y. ^ 



{i-ky. k {n-k-iy 

, (2^~l)(fc + /-2)! 

' {i-ky. 



l=k 

n-1 



l=k 



which may be proved by induction. 
Finally, with HA14|) together with 



(_2)"-i(2n-3)!! 



(A16) 



we have all ingredients of (|A1HI to calculate the a-s, the 
result being displayed in Eq. (|2SJ). The standard devia- 
tion cr„ given in (|28|l follows then from (|A12|I and ljA16p . 



APPENDIX B: LEADING PERTURBATION OF 
THE PDF FOR LARGE a FROM THE riTH MODE 

We start out from the Fourier representation H68|) of 
the path with one mode of frequency n beside the basic 
one (n = 1). From the condition h'{to) = 0, wc obtain the 
correction in the position — t:/2 + Sn of the maximum 
to leading order 



On = nSn cos — + Iflj 

ai \ 2 



(Bl) 



Hence we can calculate the maximum to second order 
in e„ (the quadratic correction in Sn contributes only to 
cubic order) 

h^ = h(tn) « h{Tr/2) + h\Tr/2)6n + l/2h"{TT/2)Sl 
« ai + £„a„ sin + 



2 

2ai \^ 



(B2) 



Now the PDF for the MRH is obtained by averaging 
over ai, a„, (p„ in 



P{z) = {Siz-h^y . 
Expanding to second order, one finds 



P{z) ^ {S{z-ai)) - {S'{z-ai) 



£„a„ sm 



/m: 



2 2 2 

„ cos"^ , „ 
2ai V 2 



+^(<5"(z-ai)a^in^(^ + ^„)) 



(B3) 



(B4) 



Note that here derivatives of the Dirac-delta appear. 
Now performing the averages yields (Pq is given by 
Eq. ^) 



P{z) = Po{z)+elP2.n{z), 



n 
Y 



0{z)ze- 



(B5) 
(B6) 



Differentiation of the terms with step-functions gives 

P2,„(2) = {l-ny2)S{z) + zS'{z)/2 

+e-''d{z) {2z^ + {n^ - 3)z) , (B7) 

where the term proportional to z'^S{z) has been omitted, 
since it does not contribute to the average and other mo- 
ments of nonsingular functions. Hence we obtain formula 

m- 

The mean to second order is best calculated from 
(|B5|) . (|B6|) resulting in 



(B8) 



The scaled PDF is then obtained by the change of vari- 
able from z to X = z/ (z), and expanding the resulting 
expression to second order in e„ yields Eqs. IjTlf) and 
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